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ABSTRACT 

The evidence for particle acceleration in supernova shells comes from electrons whose synchrotron 
emission is observed in radio and X-rays. Recent observations by the HESS instrument reveal that 
supernova remnants also emit TeV 7-rays; long awaited experimental evidence that supernova rem- 
nants can accelerate cosmic rays up to the "knee" energies. Still, uncertainty exists whether these 
7-rays are produced by electrons via inverse Compton scattering or by protons via 7r°-decay. The 
mult i- wavelength spectra of supernova remnants can be fitted with both mechanisms, although a 
preference is often given to 7r°-decay due to the spectral shape at very high energies. A recent study 
of the interstellar radiation field indicates that its energy density, especially in the inner Galaxy, is 
higher than previously thought. In this paper we evaluate the effect of the interstellar radiation field 
on the inverse Compton emission of electrons accelerated in a supernova remnant located at different 
distances from the Galactic Centre. We show that contribution of optical and infra-red photons to 
the inverse Compton emission may exceed the contribution of cosmic microwave background and in 
some cases broaden the resulted 7-ray spectrum. Additionally, we show that if a supernova remnant 
is located close to the Galactic Centre its 7-ray spectrum will exhibit a "universal" cutoff at very 
high energies due to the Klein-Nishina effect and not due to the cut-off of the electron spectrum. As 
an example, we apply our calculations to the supernova remnants RX J1713. 7-3946 and GO. 9+0.1 
recently observed by HESS. 

Subject headings: Galaxy: general — gamma-rays: theory — radiation mechanisms: non-thermal — 
ISM: cosmic rays — ISM: supernova remnants — ISM: individual (RXJ1713. 7-3946, 
GO.9+0.1) 



1. INTRODUCTION 



A new calculation (|Porter fc Stronel l2005fl of the 
Galactic interstellar radiation field (ISRF) consistent 
with multi-wavelength observations by DIRBE and FI- 
RAS indicates that the energy density of the ISRF is 
higher, particularly in the inner Galaxy, than previously 
thought. This has implications for the inverse Comp- 
ton (IC) scattering of electrons and positrons and other 
electromagnetic processes in the interstellar medium 
l|Moskalenko fc Strong 1 120001 iMoskalenkoeTall [2006). 
Another place where the enhanced ISRF may play a role 
is the IC emission off very high energy (VHE) electrons 
accelerated in a supernova remnant (SNR) shock. 

SNRs are believed to be the primary so urces of cosmic 
rays i n the Galaxy. Obser vations of X-ray ij Kc^Vjmia et^d.1 
1995) and 7-ray emission l|Aharonian et alJl2005lJ l2006) 
from SNR shocks reveal the presence of energetic par- 
ticles, thus testifying to efficient acceleration processes. 
Acceleration of particles in collisionless shocks is a mat- 
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tcr of intensive resear ch in conjunc t ion with the problem 
of co s mic ray origin llDrurvl 119831 iBlandford fc Eichlerl 
Il987t I.Tones fc Ellison! 1199 llL Current mode ls include 
nonlinear effects (e.g., Berezhko & Volk 2006) and treat 
particle acceleration using hydrodynamic codes (e.g., 
Ell ison fc Cassam-Chenail 12005]) . The predicted spec- 
trum of accelerated particles has a power-law form in 
rigidity with index which may slightly vary around —2.0. 
Conventional estimates put the maximum reachable par- 
ticle energy (protons) at or just below the knee en- 
ergy; for the case of electrons the maximum energy 
is lower due to the synchrotron and IC energy losses 
^Ellison fc Cassam-Chenail l2005). Young SNRs may be 
capable of particle acceleration up to 10 17 eV due to the 
effect of magnetic field amplification ( Bell fc LuceU2001^ 
around the shock and assuming Bohm diffusion and a Se- 
dov expansion law. 

The VHE 7-ray emission from shell-type SNRs has 
been modelled using leptonic (IC) and hadronic (n°- 
decay) scenarios (e.g., iBaring et all feOOS). The leptonic 
scenario fits the broad-band spectrum of a SNR as- 
suming a pool of accelerated electrons scattered off the 
cosmic microwave background (CMB) producing VHE 
7-rays while the magnetic field and electron spectrum 
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cut-off are tuned to fit the radio and X-ray data (e.g., 
lLazendic et~al~ll2004fl . The hadronic model fits the VHE 
7-ray spectrum assuming a beam of accelerated pro- 
tons hits a targe t, such as a nearby molecular cloud 
ijAharonianl l20A2ll . The latter, if definitively proven, 
would be the first experimental evidence of proton ac- 
celeration in SNRs. While there is no clear distinction 
between different models of VHE emission from SNRs, 
some authors tend to prefer the hadroni c scenario since it 
fits better the observed spect ral shape ijAharonian et alJ 
120061 [.Bcrezhko & Volk 12006^) . Such a preference, how- 
ever, is made based on a simplified "one-zone" model 
which typically includes CMB photons only. 

In this paper we evaluate the effect of the ISRF on the 
IC emission of electrons accelerated in SNRs located at 
different distances from the Galactic Centre (GC). As ex- 
amples, we apply our calculations to the shell-type SNR 
RX J1713.7-3946 a nd composite SNR GO.9+0.1 r ecently 
observed by HESS ijAharonian et alJl2005al I2006|) . 

2. INTERSTELLAR RADIATION FIELD 

The Galactic ISRF calculation uses a model for the 
distribution of stars in the Galaxy, a model for the dust 
distribution and properties, and a treatment of scatter- 
ing, absorption, and subsequent reemission of the stellar 
light by the dust. A brief descr i ption of our calculation 
is available in lMoskalenko et alJ l)2006|) : we re- iterate the 
main points here. 

Our stellar model assumes a type cl assification based 
on th at used in the SKY model of iWainscoat et alJ 
(1992), with modifications to account for results from 
recen t experiments such as 2MAS S, SDSS, and others 
(e.g- JOihalfeOOlt Uuric et all 120051) . and synthe tic spec- 
tral modelling studies (e.g.. iGirardi et, al]l2002l) . 

Dust is modelled with a mixture of polycyclic aromatic 
hydrocarbons, graphite, and silicate. The dust grains are 
spherical and we include details of their absorption and 
scattering efficiencies, abundances, and size distribution 
in the scattering and heating calc ulations (|Li &: Draind 
1200 It iWeingartner fc Draindl2b01|) . Stochastic and equi- 
libr ium heating of the dust grains is treated follow- 
ing iDraine fc ii l|200lD and iLi fc Drainel l|200U) . The 
dust is assumed to follow th e Galactic gas distribu- 
tion and metallic ity gradient ( Moskale nko et all 120021 
IStrong et al.ll200"H and references therein). 

Figure n shows the total ISRF at the positions we will 
calculate the IC emission in Section [3] Information on 
the calculation method, spatial and spectral distribution 
in the Galactic plane, and a brief comparison of our cal- 
culations with pre vious results and the loc ally observed 
ISRF are given in iMoskalenko et al J l)2006|) : further dis- 
cussion of the new ISRF is deferred to a forthcoming 
paper (Porter & Strong, in preparation). 

3. CALCULATIONS 

For an assumed isotropic ISRF and electron distribu- 
tion, the IC emissivity (photons cm -3 s" 1 ) is given by 

= clde 1 J d-/ e n e {e 1 )n e (j e ) '— , (1) 

where e 2 = E^/uieC 2 and e\ = e/m e c 2 are the dimen- 
sionless 7-ray and target photon energy, 7 e = E e /m e c 2 
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Fig. 1. — Interstellar radiation field energy density: solid line, 
R = kpc, z = kpc; dashed line, R = 3 kpc, z = —0.05 kpc; 
dotted line, R = 4 kpc, z = kpc; dash-dotted line, R = 7.5 kpc, 
z = kpc. 



is the electron Lorentz factor, n e and n e are the differ- 
ential target photon and elec tron number densities, and 
dQ 7 /de 2 is given by eq.(9) in l.lonesl l)1968l) . 

To illustrate the effect of different electron spectral in- 
dices on the IC emission, we take an electron number den- 
sity n e (7 e ) = tt-oTcT" 5 normalised so that J°° n e (-/ e )dj e — 

1 cm" 3 with 7 min = 1 MeV/m e c 2 for 6 = 1.8, 5 = 2.0 and 
S = 2.5. 

Figure shows the calculated IC emissivity using the 
ISRF at R — (upper) and R = 4 (lower) in the Galactic 
plane. There are two essential effects that contribute to 
produce the total emission curves in Fig. [5] First, the 
variation in the electron spectral index S — 1.8 — > 2.5 
increases the contribution by optical and infra-red (IR) 
photons to the 7-ray emission below ~ 50 GeV and ~ 1 
TeV, respectively. As S increases, the number of low en- 
ergy electrons increases; the IC scattering rate by these 
electrons off optical and IR photons becomes larger giv- 
ing more upscattered photons. In turn, as the 7-ray en- 
ergy increases, the contribution by the optical and IR 
photons to the total emissivity decreases from the reduc- 
tion of the scattering cross-section in the Klein-Nishina 
(KN) regime. Thus, softer electron spectra promote more 
efficient scattering of the optical and IR components, in- 
creasing their relative contributions. Second, the density 
of optical and IR photons varies with position in the 
Galaxy. Toward the GC, the optical photon density is 
high, almost an order of magnitude larger per target pho- 
ton energy interval when compared with the optical com- 
ponent at R = 4 kpc. Similarly, over the inner Galaxy 
the IR photon energy density is significantly higher than 
the CMB. This gives a contribution by the optical and 
IR photons to the total IC emission that is at least com- 
parable to the CMB over these regions of the Galaxy. 

Note that the marked decrease in the optical, then IR, 
then CMB emission as the 7-ray energy increases in Fig-EI 
is due to the reduction of the scattering cross-section in 
the KN regime, as described above, and not due to any 
cut-off we have introduced in the electron spectrum; if a 
cut-off was included in the electron spectrum, it would 
further decrease the contribution from the CMB relative 
to the optical and IR components (see Fig.^J). 
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Fig. 2. — Inverse Compton emissivity in the Galactic plane for 
R = kpc (top) and R = 4 kpc (bottom). Line styles: solid line, 
total; dashed, optical; dot-dashed, infra-red; dotted, CMB. For 
each panel, top-most line-sets correspond to 7.T 1 ' 8 , middle line- 
sets correspond to 7,7 2,0 and bottom-most correspond to 7,7 2 ' 5 
electron spectra respectively. 



To illustrate the effect of the ISRF in a purely leptonic 
scenario, we calculate t he synchrotron and IC emis sion 
using a one-zone model ijAharonian fc AtovanllT99^1 for 
the SNR RXJ1713.7-3946. We use an electron spectrum 
Q e (E) = Q E~ S exp (-E/E max ) where Q , 6 and E max 
are adjustable parameters. The synchrotron emission for 
a magneti c field strength B is calculated using t he usual 
formulae <|Ginzpurdll979l iGhisellini et al.lll988|) . 

Figure [S] shows the photon flux from RXJ1713. 7-3946 
calculated for a distance of d = 1 kpc (upper) and d = 6 
kpc (lower). This covers the range of distances quoted by 
vario us authors l)Lazendic e t al. 2004: Aharonia n et alJ 
2006) . The IC emission is calculated using the appropri- 
ate ISRF spectrum for each distance (see Fig. For 
both cases, we reproduce the TeV emission first: with 
S ~ 1.8-2.2 and E max ~ 15-40 TeV a reasonable 
fit can be obtained. Next, we fit the radio and X-ray 
data within the constraints imposed by the TeV emis- 
sion: S ~ 2.0 and E max ~ 15 - 25 TeV, with B ~ 5 - 20 
/iG are allowed. The electron luminosity for an age of 
1000 years (d — 1 kpc) or 10000 years (d = 6 kpc) is 
(1 — 4) x 10 37 erg s _1 with larger luminosities for the more 
distant case. The parameters used are not unique. As a 
particular example, the values used in Fig.[3]are: 6 = 2.0, 
electron luminosity of 1.5 x 10 37 erg s _1 , E max = 25 TeV, 




10' 6 10" 4 10' 2 1 10 2 10 4 10 6 10 s 10 10 10 12 10" 



E, (eV) 

Fig. 3.— Fit to the flux spectrum of RXJ1713. 7-3946 using 
the ISRF for d = 1 kpc (upper) and d = 6 kpc (lower). Line- 
styles: solid, total synchrotron and IC flux; dashed, optical IC; 
dash -dotted, infra-red IC; dotted, CMB IC. Data: radio (ATCA) 
from ILazendic e*t~ al. (2004) for the north-west ri m of RXJ1713.7- 
3946, but upscaled by a fa ctor of 2 following [Aharonian et alJ 
<2T)09): X-ray (ASCA) from IPaiiiiuti et alJ 1200311 : ^-rav (HESS1 
fromlAharonian et al. (20Q3). 



and B = 12 /j,G (d = 1 kpc); S = 2.0, electron luminosity 
of 3.5 x 10 37 erg s -1 , E max = 25 TeV, and B = 15 fiG 
(d = 6 kpc). 

We also consi der the composite SNR GO. 9+0.1 ob- 
served by HESS ij Aharonian et al.l2005bD : this composite 
SNR is a source of VHE ele ctrons, as evidenced by the 
synchrotron X-ray emission ijPorouet et a l. 2003) , and is 
located in the GC region where the ISRF density is high. 
As for RXJ1713. 7-3946, we use a one-zone model except 
we modify the source spectrum to be a broken power-law 
with a break at energy -Eb roa k- As before, we adjust the 
electron luminosity and high energy spectral index to fit 
the TeV observations, then vary B and -E max to repro- 
duce the X-ray data. Then, EWcak and the low energy 
spectral index are adjusted to give consistency with the 
radio data. Figure 01 shows the calculated intrinsic pho- 
ton flux for GO. 9+0.1 at the GC. The parameter values 
are: Sl — 2.0, 8jj — 2.8, -Ebroak = 15 GeV, where Sl 
and (5r/ are the spectral indices below and above -Ebreak 
respectively, B — 9.5 /xG, E max = 1000 TeV, and an elec- 
tron luminosity 2.5 x 10 38 erg s _1 for an assumed age of 
5000 years. Note that these values are again not unique, 
and other combinations are possible. 

The cut-off in the IC emission for the optical and IR 
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Fig. 4. — Fit to the flux spectrum of GO. 9+0.1. Line-styl es as for 
Fig. Data: radio (VLA) from Hclfand & Becker! 119871: X-ray 
(BeppoSAX) at 1.8-10 keV from lSidoli et alj 12000) (upp er solid 
region) and (XMM) at 2-12 keV fromlPorauet et^TT 120031) (lower 
solid region); 7-ray (HESS) from Aharonian et al. 1 2005a). 



in Figs. |3] and 0] comes from the reduction of the scat- 
tering cross-section in the KN regime; only for the CMB 
is the IC spectral cut-off solely due to the cut-off in the 
electron spectrum. The cut-off in the IC emission for the 
IR in Fig. is only partially due to the KN effect since 
the electron spectral cut-off energy is essentially lower. 
This effect plays an important role for a SNR close to 
the GC where the ISRF is dominant. In this case, the 
cut-off energy of the VHE 7-ray spectrum becomes uni- 
versal and defined exclusively by the reduction of the 
scattering cross-section and not by an assumed electron 
spectral cut-off. Note also, for energies higher than ~ 10 
TeV, the intrinsic spectrum will be further reduced by 
attenuation due to pair production on the ISRF itself 
(jMoskalenko et alJl2006^ . 

4. DISCUSSION 

We have evaluated the effect of the ISRF on the IC 
emission of electrons accelerated in SNRs. Our calcula- 
tion shows that for a SNR located in the inner Galaxy the 
contribution of optical and IR photons to the IC emis- 
sion exceeds the contribution of the CMB and in some 



cases broadens the resultant 7-ray spectrum. SNRs lo- 
cated close to the GC exhibit a "universal" cut-off at 
VHEs due to the KN effect and not due to a cut-off in 
the electron spectrum. 

We have made a calculation of the IC and synchrotron 
emission for simple power-law electron spectra including 
the contribution by the ISRF using the shell-type SNR 
RXJ1713.7-3946 and composite SNR GO.9+0.1 as exam- 
ples. Within the confines of a simple one-zone leptonic 
model, it is possible to fit the observed flux spectrum 
with a reasonable combination of model parameters. Ob- 
servations in the GeV to sub- TeV range by the GLAST 
experiment will be critical in distinguishing between lep- 
tonic or hadronic scenarios of 7-ray production in SNRs 
as the predictions for the spectral shape in this energy 
range are distinctly different. 

Finally, we point out at least two cases for the leptonic 
model where the contribution of the ISRF will be critical. 
First will be when fitting the synchrotron peak dictates 
a relatively low- energy cut-off in the electron spectrum. 
In this case, the VHE 7-rays will be produced by lower- 
energy electrons scattered off optical and IR photons. 
Second will be when fitting the synchrotron peak dictates 
a high- energy cut-off in the electron spectrum. An equal 
or exceeding contribution of optical and IR protons, but 
with lower energy cut-off in the spectrum of upscattered 
7-rays due to the KN effect, will effectively broaden the 
IC peak. Observations of the TeV emission from SNRs 
in the inner Galaxy may thus serve as a probe of the 
ISRF near their location. In turn, observations of shell- 
type SNRs in the outer Galaxy where the CMB photons 
will provide the majority of the IC emission may be used 
to evaluate the spectrum of accelerated electrons in the 
shell and for studies of shock acceleration. 
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